Updated 1/15/2025 by K.Lenderman - data check - added new code for 2026 data
This is an R Markdown Notebook. When you execute code within the notebook, the results appear beneath the code.
#data_all <- list.files(path = "Lab_Data_TissueProcessing/raw_data/Files_by_Month", # Identify all CSV files
# pattern = "*.csv", full.names = TRUE) %>%
#lapply(read_csv) %>% # Store all files in list
# bind_rows # Combine data sets into one data set
#data_all # Print data to RStudio console
#as.data.frame(data_all) # Convert tibble to data.frame
#Filtering NAs and unnecessary columns
#data_all <- data_all %>% filter(!is.na(date_collected))
#data_all <- select(data_all, -light_regime, -oyster_zone)
#write.csv(data_all, "Master_files/tissue_processing_all_data.csv", row.names=FALSE)
##### USE THIS CODE TO MERGE DATA FILES - some of the files have columns that are not the same format (dates specifically), which was causing issues in merging. This code below should solve that problem. It converts all date columns to dates and m/d/y format #####
file_paths <- list.files(path = "Lab_Data_TissueProcessing/raw_data/Files_by_Month",
pattern = "*.csv", full.names = TRUE)
data_all <- lapply(file_paths, function(file_path) {
read_csv(file_path) %>%
mutate(
date_collected = as.Date(date_collected,format = "%m/%d/%Y"),
date_processed = as.Date(date_processed,format = "%m/%d/%Y"),
date_davidsons = as.Date(date_davidsons,format = "%m/%d/%Y"),
date_etoh = as.Date(date_etoh,format = "%m/%d/%Y")
)
}) %>%
bind_rows
New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 38 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (8): lab_id, date_collected, site, date_processed, condition, sample_notes, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (7): mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (4): lab_id, site, oyster_zone, condition
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (20): light_regime, dissection_notes, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (6): lab_id, site, oyster_zone, light_regime, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (18): mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ...25, ...2...
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 27── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (18): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17, ...18, ...19, ...2...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 27── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (18): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17, ...18, ...19, ...2...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (5): lab_id, site, oyster_zone, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (19): light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...2...
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (6): lab_id, site, oyster_zone, light_regime, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (18): mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ...25, ...2...
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (4): lab_id, site, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (20): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22...
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (4): lab_id, site, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (20): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22...
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 27── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (18): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17, ...18, ...19, ...2...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 35 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 46 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (8): lab_id, date_collected, site, date_processed, condition, sample_notes, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (7): mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (5): lab_id, site, oyster_zone, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (19): light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...2...
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (6): lab_id, site, oyster_zone, light_regime, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (19): mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, date_etoh, ...21, ...22, ...23, ...24, ...
date (3): date_collected, date_processed, date_davidsons
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (4): lab_id, site, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (20): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22...
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (4): lab_id, site, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (20): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22...
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 27── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (8): lab_id, date_collected, site, date_processed, condition, sample_notes, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (17): mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17, ...18, ...19, ...20, ...21, ...2...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (8): lab_id, date_collected, site, date_processed, condition, sample_notes, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (7): mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (5): lab_id, site, oyster_zone, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (19): light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...2...
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 998 Columns: 33── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (6): lab_id, site, oyster_zone, light_regime, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (19): mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, date_etoh, ...21, ...22, ...23, ...24, ...
date (3): date_collected, date_processed, date_davidsons
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (4): lab_id, site, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (21): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, date_etoh, ....
date (3): date_collected, date_processed, date_davidsons
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (4): lab_id, site, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (21): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, date_etoh, ....
date (3): date_collected, date_processed, date_davidsons
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (8): lab_id, date_collected, site, date_processed, condition, sample_notes, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (7): mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (8): lab_id, date_collected, site, date_processed, condition, sample_notes, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (7): mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (8): lab_id, date_collected, site, date_processed, condition, sample_notes, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (7): mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (5): lab_id, site, oyster_zone, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (19): light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...2...
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (6): lab_id, site, oyster_zone, light_regime, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (19): mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, date_etoh, ...21, ...22, ...23, ...24, ...
date (3): date_collected, date_processed, date_davidsons
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (4): lab_id, site, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (21): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, date_etoh, ....
date (3): date_collected, date_processed, date_davidsons
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (4): lab_id, site, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (21): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, date_etoh, ....
date (3): date_collected, date_processed, date_davidsons
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (5): lab_id, site, oyster_zone, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (19): light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...2...
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (6): lab_id, site, oyster_zone, light_regime, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (18): mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ...25, ...2...
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (4): lab_id, site, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (20): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22...
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (4): lab_id, site, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (20): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22...
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (8): lab_id, date_collected, site, date_processed, condition, sample_notes, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (7): mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (5): lab_id, site, oyster_zone, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (19): light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...2...
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (6): lab_id, site, oyster_zone, light_regime, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (18): mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ...25, ...2...
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 36 Columns: 20── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (4): lab_id, site, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (7): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (4): lab_id, site, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (20): oyster_zone, light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22...
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (5): lab_id, site, oyster_zone, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (19): light_regime, mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...2...
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 33── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (6): lab_id, site, oyster_zone, light_regime, condition, dissection_notes
dbl (5): lab_sample, ww_total_g, height_mm, length_mm, width_mm
lgl (18): mantle_rftm, ggr1_etoh, ggr2_etoh, adductor_etoh, cross_histology, ...21, ...22, ...23, ...24, ...25, ...2...
date (4): date_collected, date_processed, date_davidsons, date_etoh
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 999 Columns: 27── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (8): lab_id, date_collected, site, date_processed, condition, sample_notes, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (17): mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17, ...18, ...19, ...20, ...21, ...2...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (8): lab_id, date_collected, site, date_processed, condition, sample_notes, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (7): mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.New names:Rows: 36 Columns: 17── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (8): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology, ...15, ...16, ...17
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.Rows: 9 Columns: 14── Column specification ─────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (7): lab_id, date_collected, site, date_processed, condition, date_davidsons, date_etoh
dbl (2): lab_sample, height_mm
lgl (5): sample_notes, mantle_rftm, ggr1_frozen, ggr2_etoh, cross_histology
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.
as.data.frame(data_all) # Convert tibble to data.frame
data_all <- data_all %>% filter(!is.na(date_collected))
data_all <- select(data_all, -light_regime, -oyster_zone)
data_all
NA
NA
NA
#A wrong value caliper input was identified for the height of sample 0923GOLD_23. This code is removing that value from the data set as we cannot conclude what this original value was. The value is 8.62. This will cause this individual to fall out of the dataset when standardized to length. This code does not completely remove the individual from the dataset.
data_all$height_mm[data_all$height_mm == "0"] <- NA
data_all <- data_all %>% dplyr::mutate(date_collected= as.Date(date_collected), month = month(date_collected))
data_all <- data_all %>% dplyr::mutate(date_collected= as.Date(date_collected), year = year(date_collected))
#changing numeric month to month name
data_all$month <- factor(data_all$month, levels = c("1","2","3","4","5","6","7", "8", "9", "10", "11", "12"),
labels=c("Jan","Feb", "March","April","May", "June", "July", "Aug", "Sept", "Oct", "Nov", "Dec"))
data_all
data_all<- data_all %>%
dplyr::mutate(condition_score = recode(condition), "1_very_good" = 1, "2_good" = 2, "3_good_minus"= 3, "4_fair_plus"= 4, "5_fair"= 5,"6_fair_minus"=6,"7_poor_plus"=7, "8_poor"= 8, "9_very_poor"= 9)
Error in `dplyr::mutate()`:
ℹ In argument: `condition_score = recode(condition)`.
Caused by error in `recode()`:
! argument "recodes" is missing, with no default
Run `]8;;x-r-run:rlang::last_trace()rlang::last_trace()]8;;` to see where the error occurred.
#This chunk of code is removing 0723LAUR_20 and 0723LAUR_26 from the datasheet as they have been identified as spat on shell to avoid bias in the data. During this sample collection there were animals that were significantly smaller than the single set oysters. These individuals should be removed from all monthly sampling related datasheets including disease analysis. All tissue amples will be disgarded.
data_all <- data_all %>%
subset(lab_id != "0723LAUR_20") %>%
subset(lab_id != "0723LAUR_26")
data_all
NA
st_height <- summarySE(data_all%>% filter(!is.na(height_mm)), measurevar="height_mm", groupvars=c("site", "date_collected"))
st_height
#Calculate completeness for QC
st_height$Completeness <- st_height$N /30
st_height
write.csv(st_height, "Lab_Data_TissueProcessing/output\\Completeness_tissue_processing_data.csv", row.names=FALSE)
#Mean Height
ggplot(data=data_all, aes(x=site, y=height_mm, fill=site)) +
geom_boxplot()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title="Mean Shell Height ", x ="site", y = "Mean Shell Height (mm)") + facet_wrap(.~year)
#Body condition
ggplot(data=data_all, aes(x=site, y=condition_score, fill=site)) +
geom_boxplot()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title="Mean Condition Score ", x ="site", y = "Mean Body Condition Score")+ scale_y_reverse()+ facet_wrap(.~year)
mean_body_condition
Error: object 'mean_body_condition' not found
#Proportions graph Body condition scores
#Proportions graph Body condition scores - 2023
df_BCS_proportions2023<- data_2023 %>%
mutate(condition_score_numeric = as.numeric(condition_score),condition_score_bin = case_when(condition_score == 1 ~ "1", condition_score == 2 ~"2",condition_score == 3 ~ "3", condition_score == 4 ~"4",condition_score == 5 ~ "5", condition_score == 6 ~"6",condition_score == 7 ~ "7", condition_score == 8 ~"8", condition_score ==9 ~"9", TRUE ~ as.character(condition_score))) %>%
group_by(site, month, year, condition_score_bin) %>% dplyr ::summarise(Count= n()) %>%
ungroup() %>%
mutate(Proportion = Count/sum(Count))
`summarise()` has grouped output by 'site', 'month', 'year'. You can override using the `.groups` argument.
df_BCS_proportions2023<- na.omit(df_BCS_proportions2023)
BCS_proportion_all_2023<- ggplot(data=df_BCS_proportions2023, aes(x=month, y= Proportion, fill=condition_score_bin)) +
geom_bar(width = .5, stat="identity", position = "fill", colour = "black")+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1))+
labs(title="Proportion of body condition scores - 2023", x ="month", y = "Proportion of body condition scores ")+ theme(axis.title.y = element_text(size = rel(1.3), angle =90), axis.title.x = element_text(size = rel(1.3), angle = 0))+
theme(axis.text=element_text(size=12))+
scale_fill_brewer(palette = "Blues", direction = -1)+
facet_wrap(~ site)
BCS_proportion_all_2023
#pdf(paste0(path = "Lab_Data_TissueProcessing/output" ,"/all_BCS_proportion_2023.pdf"), height = 7, width = 13)
#print(BCS_proportion_all_2023)
#dev.off()
#Proportions graph Body condition scores - 2024
df_BCS_proportions2024<- data_2024 %>%
mutate(condition_score_numeric = as.numeric(condition_score),condition_score_bin = case_when(condition_score == 1 ~ "1", condition_score == 2 ~"2",condition_score == 3 ~ "3", condition_score == 4 ~"4",condition_score == 5 ~ "5", condition_score == 6 ~"6",condition_score == 7 ~ "7", condition_score == 8 ~"8", condition_score ==9 ~"9", TRUE ~ as.character(condition_score))) %>%
group_by(site, month, year, condition_score_bin) %>% dplyr ::summarise(Count= n()) %>%
ungroup() %>%
mutate(Proportion = Count/sum(Count))
`summarise()` has grouped output by 'site', 'month', 'year'. You can override using the `.groups` argument.
df_BCS_proportions2024<- na.omit(df_BCS_proportions2024)
BCS_proportion_all_2024<- ggplot(data=df_BCS_proportions2024, aes(x=month, y= Proportion, fill=condition_score_bin)) +
geom_bar(width = .5, stat="identity", position = "fill", colour = "black")+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1))+
labs(title="Proportion of body condition scores - 2024", x ="month", y = "Proportion of body condition scores ")+ theme(axis.title.y = element_text(size = rel(1.3), angle =90), axis.title.x = element_text(size = rel(1.3), angle = 0))+
theme(axis.text=element_text(size=12))+
scale_fill_brewer(palette = "Blues", direction = -1)+
facet_wrap(~ site)
BCS_proportion_all_2024
#pdf(paste0(path = "Lab_Data_TissueProcessing/output" ,"/all_BCS_proportion_2024.pdf"), height = 7, width = 13)
#print(BCS_proportion_all_2024)
#dev.off()
#Proportions graph Body condition scores - 2025
df_BCS_proportions2025<- data_2025 %>%
mutate(condition_score_numeric = as.numeric(condition_score),condition_score_bin = case_when(condition_score == 1 ~ "1", condition_score == 2 ~"2",condition_score == 3 ~ "3", condition_score == 4 ~"4",condition_score == 5 ~ "5", condition_score == 6 ~"6",condition_score == 7 ~ "7", condition_score == 8 ~"8", condition_score ==9 ~"9", TRUE ~ as.character(condition_score))) %>%
group_by(site, month, year, condition_score_bin) %>% dplyr ::summarise(Count= n()) %>%
ungroup() %>%
mutate(Proportion = Count/sum(Count))
`summarise()` has grouped output by 'site', 'month', 'year'. You can override using the `.groups` argument.
df_BCS_proportions2025<- na.omit(df_BCS_proportions2025)
BCS_proportion_all_2025<- ggplot(data=df_BCS_proportions2025, aes(x=month, y= Proportion, fill=condition_score_bin)) +
geom_bar(width = .5, stat="identity", position = "fill", colour = "black")+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1))+
labs(title="Proportion of body condition scores - 2025", x ="month", y = "Proportion of body condition scores ")+ theme(axis.title.y = element_text(size = rel(1.3), angle =90), axis.title.x = element_text(size = rel(1.3), angle = 0))+
theme(axis.text=element_text(size=12))+
scale_fill_brewer(palette = "Blues", direction = -1)+
facet_wrap(~ site)#, scales = "free")
BCS_proportion_all_2025
# pdf(paste0(path = "Lab_Data_TissueProcessing/output" ,"/all_BCS_proportion_2025.pdf"), height = 7, width = 13)
# print(BCS_proportion_all_2025)
# dev.off()
#Proportions graph Body condition scores - 2026
df_BCS_proportions2026<- data_2026 %>%
mutate(condition_score_numeric = as.numeric(condition_score),condition_score_bin = case_when(condition_score == 1 ~ "1", condition_score == 2 ~"2",condition_score == 3 ~ "3", condition_score == 4 ~"4",condition_score == 5 ~ "5", condition_score == 6 ~"6",condition_score == 7 ~ "7", condition_score == 8 ~"8", condition_score ==9 ~"9", TRUE ~ as.character(condition_score))) %>%
group_by(site, month, year, condition_score_bin) %>% dplyr ::summarise(Count= n()) %>%
ungroup() %>%
mutate(Proportion = Count/sum(Count))
`summarise()` has grouped output by 'site', 'month', 'year'. You can override using the `.groups` argument.
df_BCS_proportions2026<- na.omit(df_BCS_proportions2026)
BCS_proportion_all_2026<- ggplot(data=df_BCS_proportions2026, aes(x=month, y= Proportion, fill=condition_score_bin)) +
geom_bar(width = .5, stat="identity", position = "fill", colour = "black")+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1))+
labs(title="Proportion of body condition scores - 2026", x ="Month", y = "Proportion of body condition scores ")+ theme(axis.title.y = element_text(size = rel(1.3), angle =90), axis.title.x = element_text(size = rel(1.3), angle = 0))+
theme(axis.text=element_text(size=12))+
scale_fill_brewer(palette = "Blues", direction = -1)+
facet_wrap(~ site)#, scales = "free")
BCS_proportion_all_2026
pdf(paste0(path = "Lab_Data_TissueProcessing/output" ,"/all_BCS_proportion_2026.pdf"), height = 7, width = 13)
print(BCS_proportion_all_2026)
dev.off()
png
2
#Winter sampling
# winter_bcs <- rbind(df_BCS_proportions2024, df_BCS_proportions2025)
#
# winter_bcs <- winter_bcs %>% filter(month == "Dec"|month == "Jan"|month == "Feb")
#
# BCS_proportion_winter<- ggplot(data=winter_bcs, aes(x=factor (month, level=c('Dec', 'Jan', 'Feb')), y= Proportion, fill=condition_score_bin)) +
# geom_bar(width = .5, stat="identity", position = "fill", colour = "black")+
# theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
# theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1))+
# labs(title="Proportion of body condition scores - Winter Sampling 2024-2025", x ="Month", y = "Proportion of body condition scores ")+ theme(axis.title.y = element_text(size = rel(1.3), angle =90), axis.title.x = element_text(size = rel(1.3), angle = 0))+
# theme(axis.text=element_text(size=12))+
# scale_fill_brewer(palette = "Blues", direction = -1)+
# facet_wrap(~ site)
#
# BCS_proportion_winter
# Proportion per month of intensity scores at three sites during 2023 & 2024 & 2025
df_bcs_proportions_new <- rbind(df_BCS_proportions2023, df_BCS_proportions2024, df_BCS_proportions2025, df_BCS_proportions2026)
df_bcs_proportions_new<- df_bcs_proportions_new %>% filter(!site=="LAUR")
ggplot(data=df_bcs_proportions_new, aes(x=month, y= Proportion, fill=condition_score_bin)) +
geom_bar(width = .5, stat="identity", position = "fill",colour = "black")+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1))+
labs(fill ="Body condition scores")+
labs(title="Proportion of body condition scores", x ="month", y = "Proportion")+ theme(axis.title.y = element_text(size = rel(1.3), angle =90), axis.title.x = element_text(size = rel(1.3), angle = 0))+
theme(axis.text=element_text(size=12))+
scale_fill_brewer(palette = "Blues", direction = -1)+
#scale_x_continuous("Month", breaks = c(3,4,5,6,7,8,9,10))+
facet_grid(year ~site)
#Shell Pathology - all data #Includes 2023 data and all data that is in the 2024 format
shell_path <- list.files(path = "Lab_Data_TissueProcessing/raw_data/shell_pathology", # Identify all CSV files
pattern = "*.csv", full.names = TRUE) %>%
lapply(read_csv) %>% # Store all files in list
bind_rows # Combine data sets into one data set
shell_path <- shell_path %>% filter(!is.na(lab_id))
shell_path <- shell_path %>% select(-c("...14":"...26"))
write.csv(shell_path, "Master_files/shell_pathology_all_data.csv", row.names=FALSE)
#Writing summary for shell pathology
# Clean up, convert to integer, drop ID, etc.
shell_path_summary <- shell_path %>%
# split `lab_id` into two parts
separate(lab_id, into = c("Date_Site", "ID")) %>%
# drop unwanted columns
select(-pathology_notes, -ID) %>%
# convert a set of logical (or whatever) columns to 0/1
mutate(across(
c(boring_sponge, polydora, conchiolin_mod_severe, mud_blister,
pea_crab, gill_erosion, pale_digestive, discoloration,
horn_add, cyst_abscess, tumor, oyster_drill,
boring_snail, shell_scarring),
~ as.integer(as.logical(.))
))
# Compute a pathogen‐indicator sum *per sample* (row).
shell_path_summary <- shell_path_summary %>%
mutate(
Pathogen_count = rowSums(
select(., boring_sponge, polydora, conchiolin_mod_severe, mud_blister,
pea_crab, gill_erosion, pale_digestive, discoloration,
horn_add, cyst_abscess, tumor, oyster_drill,
boring_snail, shell_scarring),
na.rm = TRUE
)
)
# Now group + summarise
shell_path_count <- shell_path_summary %>%
group_by(Date_Site) %>%
summarise(
Sample_count = n(),
across(
c(boring_sponge, polydora, conchiolin_mod_severe, mud_blister,
pea_crab, gill_erosion, pale_digestive, discoloration,
horn_add, cyst_abscess, tumor, oyster_drill,
boring_snail, shell_scarring, Pathogen_count),
sum,
na.rm = TRUE
)
) %>%
ungroup()
shell_path_count
shell_path_summary
write.csv(shell_path_count, "Lab_Data_TissueProcessing/output/shell_pathology_counts.csv", row.names=FALSE)
Individual site summaries and QAQC
df_ASHC<- data_all%>%
filter(site=="ASHC")
df_ASHC
## Shell Height ##
st_height_ASHC <- summarySE(df_ASHC, measurevar="height_mm", groupvars=c("date_collected"))
st_height_ASHC
## Body Condition ##
#Excludes April and May due to scoring change. These months are scored categorically 'fat, medium, watery'.
st_condition_ASHC<- summarySE(df_ASHC, measurevar = "condition_score", groupvars = c("date_collected"))
st_condition_ASHC
mean_condition_ASHC <- st_condition_ASHC %>%
filter(!is.na(condition_score)) %>% # Filter out rows where condition_score is NA
group_by(date_collected)%>%
summarize(mean_bcs = mean(condition_score, na.rm = TRUE))
mean_condition_ASHC
ggplot(data=df_ASHC, aes(x=month, y=height_mm, group= month)) +
geom_boxplot()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title=" Ash Creek Mean Shell Height ", x ="month", y = "Mean Shell Height (mm)")+facet_wrap(.~year)
#Condition distribution across all sample months - 2023
df_ASHC_2023 <- df_ASHC%>% filter(!year=="2024")
ASHC_BCS_dist2023 <-ggplot(data=df_ASHC_2023, aes(x= condition)) +
geom_bar()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title="Ash Creek Body Condition Score - 2023", x ="Condition categorization")+
facet_wrap(~ month, scales = "free")
ASHC_BCS_dist2023
#pdf(paste0(path = "Lab_Data_TissueProcessing/output" ,"/ASHC_BCS_dist_2023.pdf"), height = 7, width = 13)
#print(ASHC_BCS_dist2023)
#dev.off()
#Condition distribution across all sample months - 2024
df_ASHC_2024 <- df_ASHC%>% filter(!year=="2023")
ASHC_BCS_dist2024 <-ggplot(data=df_ASHC_2024, aes(x= condition)) +
geom_bar()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title="Ash Creek Body Condition Score - 2024", x ="Condition categorization")+
facet_wrap(~ month, scales = "free")
ASHC_BCS_dist2024
#pdf(paste0(path = "Lab_Data_TissueProcessing/output" ,"/ASHC_BCS_dist_2024.pdf"), height = 7, width = 13)
#print(ASHC_BCS_dist2024)
#dev.off()
#Condition distribution across all sample months - 2025
df_ASHC_2025 <- df_ASHC%>% filter(year=="2025")
ASHC_BCS_dist2025 <-ggplot(data=df_ASHC_2025, aes(x= condition)) +
geom_bar()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title="Ash Creek Body Condition Score - 2025", x ="Condition categorization")+
facet_wrap(~ month, scales = "free")
ASHC_BCS_dist2025
# pdf(paste0(path = "Lab_Data_TissueProcessing/output" ,"/ASHC_BCS_dist_2025.pdf"), height = 7, width = 13)
# print(ASHC_BCS_dist2025)
# dev.off()
#Condition distribution across all sample months - 2026
df_ASHC_2026 <- df_ASHC%>% filter(year=="2026")
ASHC_BCS_dist2026 <-ggplot(data=df_ASHC_2026, aes(x= condition)) +
geom_bar()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title="Ash Creek Body Condition Score - 2026", x ="Condition categorization")+
facet_wrap(~ month, scales = "free")
ASHC_BCS_dist2026
pdf(paste0(path = "Lab_Data_TissueProcessing/output" ,"/ASHC_BCS_dist_2026.pdf"), height = 7, width = 13)
print(ASHC_BCS_dist2026)
dev.off()
#Mean Body condition per month
ggplot(data=df_ASHC, aes(x= month, y= condition_score, group = month)) +
geom_boxplot()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title="Ash Creek Mean Body Condition Score ", x ="month", y= " condition score (1-9)") + scale_y_reverse()+facet_wrap(.~year)
#ASHC Proportions graph Body condition scores - all years
ASHC_BCS_proportions<- df_ASHC %>%
mutate(condition_score_numeric = as.numeric(condition_score),condition_score_bin = case_when(condition_score == 1 ~ "1", condition_score == 2 ~"2",condition_score == 3 ~ "3", condition_score == 4 ~"4",condition_score == 5 ~ "5", condition_score == 6 ~"6",condition_score == 7 ~ "7", condition_score == 8 ~"8", condition_score ==9 ~"9", TRUE ~ as.character(condition_score))) %>%
group_by(site, month, condition_score_bin, year) %>% dplyr::summarise(Count= n()) %>%
ungroup() %>%
mutate(Proportion = Count/sum(Count))
ASHC_BCS_proportions<- na.omit(ASHC_BCS_proportions)
ASHC_BCS_proportions
BCS_proportion_ASHC<- ggplot(data=ASHC_BCS_proportions, aes(x=month, y= Proportion, fill=condition_score_bin)) +
geom_bar(width = .5, stat="identity", position = "fill")+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1))+
labs(title="Ash Creek Proportion of body condition scores June- November", x ="month", y = "Proportion of body condition scores ")+ theme(axis.title.y = element_text(size = rel(1.3), angle =90), axis.title.x = element_text(size = rel(1.3), angle = 0))+
theme(axis.text=element_text(size=12))+
facet_wrap(~ year)
#scale_fill_brewer() +
#facet_wrap(~ site)
BCS_proportion_ASHC
pdf(paste0(path = "Lab_Data_TissueProcessing/output" ,"/ASHC_BCS_proportion.pdf"), height = 7, width = 13)
print(BCS_proportion_ASHC)
dev.off()
df_ASHC
ASHC_percent_greater_3 <- df_ASHC %>%
dplyr::group_by(month, site, year) %>%
dplyr::summarise(Percentage = mean(condition_score <= 3)*100)
ASHC_percent_greater_3
ASHC_BCS_percentage <- ASHC_percent_greater_3%>%
#filter(year =="2024")%>%
ggplot(aes(x = month, y = Percentage)) +
#geom_bar(width = 0.5, stat = "identity", position = "fill") +
geom_col()+
theme_bw() +
theme(
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
axis.text.x = element_text(angle = 45, vjust = 1, hjust = 1),
axis.title.y = element_text(size = rel(1.3), angle = 90),
axis.title.x = element_text(size = rel(1.3), angle = 0),
axis.text = element_text(size = 12)
) +
labs(
title = "Ash Creek % of body condition scores >= 3",
x = "month",
y = "Percentage of body condition scores >= 3"
) +
facet_wrap(~year)
#scale_x_continuous( breaks = seq(5,12, by =1) )+
#ylim(0,60)+
#scale_y_continuous(limits = c(0,100), breaks = seq(0,100, by = 10))
ASHC_BCS_percentage
pdf(paste0(path = "Lab_Data_TissueProcessing/output", "/ASHC_BCS_percentage.pdf"),height = 7, width = 13)
print(ASHC_BCS_percentage)
dev.off()
df_FENC<- data_all%>%
filter(site=="FENC")
df_FENC
## Shell Height ##
st_height_FENC <- summarySE(df_FENC, measurevar="height_mm", groupvars=c("date_collected"))
st_height_FENC
## Body Condition ##
#Excludes April and May due to scoring change. These months are scored categorically 'fat, medium, watery'.
st_condition_FENC<- summarySE(df_FENC%>% filter(!is.na(condition_score)), measurevar = "condition_score", groupvars = c("date_collected"))
st_condition_FENC
ggplot(data=df_FENC, aes(x=month, y=height_mm, group=month)) +
geom_boxplot()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title=" Fence Creek Mean Shell Height ", x ="month", y = "Mean Shell Height (mm)")+ facet_wrap(.~year)
#Condition distribution across all sample months - 2023
df_FENC.2023 <- df_FENC%>% filter(!year=="2024")
ggplot(data=df_FENC.2023, aes(x= condition)) +
geom_bar()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title="Fence Creek Body Condition index - 2023", x ="Condition categorization")+
facet_wrap(~ month, scales = "free")
#Condition distribution across all sample months - 2024
df_FENC.2024 <- df_FENC%>% filter(year=="2024")
ggplot(data=df_FENC.2024, aes(x= condition)) +
geom_bar()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title="Fence Creek Body Condition index - 2024", x ="Condition categorization")+
facet_wrap(~ month, scales = "free")
#Condition distribution across all sample months - 2025
df_FENC.2025 <- df_FENC%>% filter(year=="2025")
ggplot(data=df_FENC.2025, aes(x= condition)) +
geom_bar()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title="Fence Creek Body Condition index - 2025", x ="Condition categorization")+
facet_wrap(~ month, scales = "free")
#Condition distribution across all sample months - 2026
df_FENC.2026 <- df_FENC%>% filter(year=="2026")
ggplot(data=df_FENC.2026, aes(x= condition)) +
geom_bar()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title="Fence Creek Body Condition index - 2026", x ="Condition categorization")+
facet_wrap(~ month, scales = "free")
#Mean Body condition per month
ggplot(data=df_FENC, aes(x= month, y= condition_score, group = month)) +
geom_boxplot()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title="Fence Creek Mean Body Condition Score", x ="month", y= " condition score (1-9)")+ scale_y_reverse() + facet_wrap(~ year, scales = "free")
#FENC Proportions graph Body condition scores - all years
FENC_BCS_proportions<- df_FENC %>%
mutate(condition_score_numeric = as.numeric(condition_score),condition_score_bin = case_when(condition_score == 1 ~ "1", condition_score == 2 ~"2",condition_score == 3 ~ "3", condition_score == 4 ~"4",condition_score == 5 ~ "5", condition_score == 6 ~"6",condition_score == 7 ~ "7", condition_score == 8 ~"8", condition_score ==9 ~"9", TRUE ~ as.character(condition_score))) %>%
group_by(site, month, condition_score_bin, year) %>% dplyr::summarise(Count= n()) %>%
ungroup() %>%
mutate(Proportion = Count/sum(Count))
FENC_BCS_proportions<- na.omit(FENC_BCS_proportions)
FENC_BCS_proportions
BCS_proportion_FENC<- ggplot(data=FENC_BCS_proportions, aes(x=month, y= Proportion, fill=condition_score_bin)) +
geom_bar(width = .5, stat="identity", position = "fill")+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1))+
labs(title="Fence Creek Proportion of body condition scores", x ="month", y = "Proportion of body condition scores ")+ theme(axis.title.y = element_text(size = rel(1.3), angle =90), axis.title.x = element_text(size = rel(1.3), angle = 0))+
theme(axis.text=element_text(size=12))+
facet_wrap(~ year)
#scale_fill_brewer() +
#facet_wrap(~ site)
BCS_proportion_FENC
pdf(paste0(path = "Lab_Data_TissueProcessing/output" ,"/FENC_BCS_proportion.pdf"), height = 7, width = 13)
print(BCS_proportion_FENC)
dev.off()
#removing 0823FENC_28 bcs was NA
df_FENC <- df_FENC %>% drop_na(condition)
FENC_percent_greater_3 <- df_FENC %>%
dplyr::group_by(month, year) %>%
dplyr::summarise(Percentage = mean(condition_score <= 3)*100)
FENC_percent_greater_3
FENC_BCS_percentage <- FENC_percent_greater_3 %>%
filter(year =="2025")%>%
ggplot(aes(x = month, y = Percentage)) +
#geom_bar(width = 0.5, stat = "identity", position = "fill") +
geom_col()+
theme_bw() +
theme(
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
axis.text.x = element_text(angle = 45, vjust = 1, hjust = 1),
axis.title.y = element_text(size = rel(1.3), angle = 90),
axis.title.x = element_text(size = rel(1.3), angle = 0),
axis.text = element_text(size = 12)
) +
labs(
title = "Fence Creek % of body condition scores >= 3",
x = "month",
y = "Percentage of body condition scores >= 3"
) +
facet_wrap(~year)
#scale_x_continuous( breaks = seq(5,12, by =1) )+
#ylim(0,60)+
#scale_y_continuous(limits = c(0,100), breaks = seq(0,100, by = 10))
FENC_BCS_percentage
pdf(paste0(path = "Lab_Data_TissueProcessing/output", "/FENC_BCS_percentage.pdf"),height = 7, width = 13)
print(FENC_BCS_percentage)
dev.off()
df_GOLD<- data_all%>%
filter(site=="GOLD")
df_GOLD <- df_GOLD %>% filter(!is.na(height_mm))
#2023 data
df_GOLD2023 <- df_GOLD%>% filter(!year=="2024")
#2024 data
df_GOLD2024 <- df_GOLD%>% filter(!year=="2023")
df_GOLD$height_mm<-as.numeric(df_GOLD$height_mm)
mean_shell_height <- df_GOLD %>% mutate(year = year(date_collected)) %>% group_by(year)%>% dplyr::summarise(mean_height = mean(height_mm)) #summarySE(measurevar="height_mm", groupvars=c("year"))
mean_shell_height
## Shell Height ##
st_height_GOLD <- summarySE(df_GOLD, measurevar="height_mm", groupvars=c("date_collected"))
st_height_GOLD
## Body Condition ##
st_condition_GOLD<- summarySE(df_GOLD, measurevar = "condition_score", groupvars = c("date_collected"))
st_condition_GOLD
#GOLD Proportions graph Body condition scores
GOLD_BCS_proportions<- df_GOLD %>%
mutate(condition_score_numeric = as.numeric(condition_score),condition_score_bin = case_when(condition_score == 1 ~ "1", condition_score == 2 ~"2",condition_score == 3 ~ "3", condition_score == 4 ~"4",condition_score == 5 ~ "5", condition_score == 6 ~"6",condition_score == 7 ~ "7", condition_score == 8 ~"8", condition_score ==9 ~"9", TRUE ~ as.character(condition_score))) %>%
group_by(site, month, year, condition_score_bin) %>% dplyr::summarise(Count= n()) %>%
ungroup() %>%
mutate(Proportion = Count/sum(Count))
GOLD_BCS_proportions<- na.omit(GOLD_BCS_proportions)
BCS_proportion_GOLD <- ggplot(data=GOLD_BCS_proportions, aes(x=month, y= Proportion, fill=condition_score_bin)) +
geom_bar(width = .5, stat="identity", position = "fill", colour = "black")+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1))+
labs(title="Gold Star Beach Proportion of body condition scores", x ="month", y = "Proportion of body condition scores ")+ theme(axis.title.y = element_text(size = rel(1.3), angle =90), axis.title.x = element_text(size = rel(1.3), angle = 0))+
theme(axis.text=element_text(size=12))+
scale_fill_brewer(palette = "Blues", direction = -1)+
facet_wrap(~ year)
BCS_proportion_GOLD
pdf(paste0(path = "Lab_Data_TissueProcessing/output" ,"/GOLD_BCS_proportion.pdf"), height = 7, width = 13)
print(BCS_proportion_GOLD)
dev.off()
df_GOLD
GOLD_percent_greater_3 <- df_GOLD %>%
dplyr::group_by(month, year) %>%
dplyr::summarise(Percentage = mean(condition_score <= 3)*100)
GOLD_percent_greater_3
GOLD_BCS_percentage <- GOLD_percent_greater_3 %>%
ggplot(aes(x = month, y = Percentage)) +
#geom_bar(width = 0.5, stat = "identity", position = "fill") +
geom_col()+
theme_bw() +
theme(
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
axis.text.x = element_text(angle = 45, vjust = 1, hjust = 1),
axis.title.y = element_text(size = rel(1.3), angle = 90),
axis.title.x = element_text(size = rel(1.3), angle = 0),
axis.text = element_text(size = 12)
) +
labs(
title = "Gold Star Beach % of body condition scores >= 3",
x = "month",
y = "Percentage of body condition scores >= 3"
) +facet_wrap(~year) +
#scale_x_continuous( breaks = seq(5,12, by =1) )+
#ylim(0,100)+
scale_y_continuous(limits = c(0,100), breaks = seq(0,100, by = 10))
GOLD_BCS_percentage
pdf(paste0(path = "Lab_Data_TissueProcessing/output", "/GOLD_BCS_percentage.pdf"),height = 7, width = 13)
print(GOLD_BCS_percentage)
dev.off()
ggplot(data=df_GOLD, aes(x=month, y=height_mm, group= month)) +
geom_boxplot()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title="Goldstar Beach Mean Shell Height ", x ="month", y = "Mean Shell Height (mm)") + #scale_x_continuous(limits= c(4,11), breaks = seq(5,10, by =1))+
geom_smooth(method = "lm", se = FALSE)+ facet_wrap(~year)
#Condition distribution across all sample months - 2023
ggplot(data=df_GOLD2023, aes(x= condition)) +
geom_bar()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title="Goldstar Beach Body Condition Score - 2023 ", x ="Condition categorization")+
facet_wrap(~ month, scales = "free")
#Condition distribution across all sample months - 2024
ggplot(data=df_GOLD2024, aes(x= condition)) +
geom_bar()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title="Goldstar Beach Body Condition Score - 2024", x ="Condition categorization")+
facet_wrap(~ month, scales = "free")
#Mean Body condition per month
ggplot(data=df_GOLD, aes(x= month, y= condition_score, group = month)) +
geom_boxplot()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
labs(title="Goldstar Beach Mean Body Condition Score ", x ="month", y= " condition score (1-9)")+ scale_y_reverse() + facet_wrap(~year)
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
df_LAUR<- data_all%>%
filter(site=="LAUR")
## Shell Height ##
st_height_LAUR <- summarySE(df_LAUR, measurevar="height_mm", groupvars=c("date_collected"))
st_height_LAUR
## Body Condition ##
## Excludes May due to scoring system change. May scored fat, medium, watery.
st_condition_LAUR<- summarySE(df_LAUR, measurevar = "condition_score", groupvars = c("date_collected"))
st_condition_LAUR
ggplot(data=df_LAUR, aes(x=month, y=height_mm, group= month)) +
geom_boxplot()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title=" Laurel Hollow Mean Shell Height - 2023", x ="month", y = "Mean Shell Height (mm)")
#Condition distribution across all sample months
ggplot(data=df_LAUR, aes(x= condition)) +
geom_bar()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
labs(title="Laurel Hollow Body Condition index - 2023", x ="Condition categorization")+
facet_wrap(~ month, scales = "free")
#Mean Body condition per month
ggplot(data=df_LAUR, aes(x= month, y= condition_score, group = month)) +
geom_boxplot()+ #scale_fill_manual(values=c("forestgreen","orange", "purple"))+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
labs(title="Laurel Hollow Mean Body Condition Score - 2023", x ="month", y= " condition score (1-9)")+ scale_y_reverse()
#theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))+
#LAUR Proportions graph Body condition scores - 2023 (only sampled one year)
LAUR_BCS_proportions<- df_LAUR %>%
mutate(condition_score_numeric = as.numeric(condition_score),condition_score_bin = case_when(condition_score == 1 ~ "1", condition_score == 2 ~"2",condition_score == 3 ~ "3", condition_score == 4 ~"4",condition_score == 5 ~ "5", condition_score == 6 ~"6",condition_score == 7 ~ "7", condition_score == 8 ~"8", condition_score ==9 ~"9", TRUE ~ as.character(condition_score))) %>%
group_by(site, month, condition_score_bin, year) %>% dplyr::summarise(Count= n()) %>%
ungroup() %>%
mutate(Proportion = Count/sum(Count))
LAUR_BCS_proportions<- na.omit(LAUR_BCS_proportions)
LAUR_BCS_proportions
BCS_proportion_LAUR<- ggplot(data=LAUR_BCS_proportions, aes(x=month, y= Proportion, fill=condition_score_bin)) +
geom_bar(width = .5, stat="identity", position = "fill")+
theme_bw() + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1))+
labs(title="Laurel Hollow Proportion of body condition scores - 2023", x ="month", y = "Proportion of body condition scores ")+ theme(axis.title.y = element_text(size = rel(1.3), angle =90), axis.title.x = element_text(size = rel(1.3), angle = 0))+
theme(axis.text=element_text(size=12))
#facet_wrap(~ year)
#scale_fill_brewer() +
#facet_wrap(~ site)
BCS_proportion_LAUR
pdf(paste0(path = "Lab_Data_TissueProcessing/output" ,"/LAUR_BCS_proportion.pdf"), height = 7, width = 13)
print(BCS_proportion_LAUR)
dev.off()
#removing 0823FENC_28 bcs was NA
df_LAUR. <- df_LAUR%>%filter(!row_number() %in% c(2))
LAUR_percent_greater_3 <- df_LAUR. %>%
dplyr::group_by(month, site, year) %>%
dplyr::summarise(Percentage = mean(condition_score <= 3)*100)
LAUR_percent_greater_3
LAUR_BCS_percentage <- LAUR_percent_greater_3 %>%
ggplot(aes(x = month, y = Percentage)) +
#geom_bar(width = 0.5, stat = "identity", position = "fill") +
geom_col()+
theme_bw() +
theme(
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
axis.text.x = element_text(angle = 45, vjust = 1, hjust = 1),
axis.title.y = element_text(size = rel(1.3), angle = 90),
axis.title.x = element_text(size = rel(1.3), angle = 0),
axis.text = element_text(size = 12)
) +
labs(
title = "Gold Star Beach % of body condition scores >= 3",
x = "month",
y = "Percentage of body condition scores >= 3"
) #+facet_wrap(~year)
#+ scale_x_continuous( breaks = seq(5,12, by =1) )+
#ylim(0,60)+
#scale_y_continuous(limits = c(0,100), breaks = seq(0,100, by = 10))
LAUR_BCS_percentage
pdf(paste0(path = "Lab_Data_TissueProcessing/output", "/LAUR_BCS_percentage.pdf"),height = 7, width = 13)
print(LAUR_BCS_percentage)
dev.off()
#Coding graphs for LIS Conference - only 2023 data
data_2023 <- data_all%>% filter(!year=="2024")
df_BCS_proportions<- data_2023 %>%
mutate(condition_score_numeric = as.numeric(condition_score),condition_score_bin = case_when(condition_score == 1 ~ "1", condition_score == 2 ~"2",condition_score == 3 ~ "3", condition_score == 4 ~"4",condition_score == 5 ~ "5", condition_score == 6 ~"6",condition_score == 7 ~ "7", condition_score == 8 ~"8", condition_score ==9 ~"9", TRUE ~ as.character(condition_score))) %>%
group_by(site, month, condition_score_bin) %>% dplyr ::summarise(Count= n()) %>%
ungroup() %>%
mutate(Proportion = Count/sum(Count))
df_BCS_proportions<- na.omit(df_BCS_proportions)
df_BCS_proportions$month <- factor(df_BCS_proportions$month, levels = c("5","6","7", "8", "9", "10", "11"),
labels=c("May", "June", "July", "Aug", "Sept", "Oct", "Nov"))
df_BCS_proportions$site <- factor(df_BCS_proportions$site, levels = c("ASHC","FENC","GOLD", "LAUR"),
labels=c("Ash Creek", "Fence Creek", "Gold Star", "Laurel Hollow"))
bcs.prop.2023 <- ggplot(data=df_BCS_proportions, aes(x=month, y= Proportion, fill=condition_score_bin)) +
geom_bar(width = .5, stat="identity", position = "fill", colour = "black")+
theme_bw() + theme(panel.grid.minor = element_blank())+
theme(axis.text.x = element_text(angle = 45, vjust = 1, hjust=1, size = rel(.9)))+
theme(panel.grid.major.x =element_blank())+
labs(title="Body condition scores in 2023", x ="Month", y = "Proportion of body condition scores ")+ theme(axis.title.y = element_text(size = rel(1.2), angle =90), axis.title.x = element_text(size = rel(1.2), angle = 0))+ scale_fill_brewer(palette = "Blues") +
theme(strip.text = element_text(size = 15))+
theme (legend.title = element_text(size = 18))+
theme(title = element_text(size = 17))+
theme(axis.text=element_text(size=17))+
labs(fill ="Body Condition Score")+
facet_wrap(~ site )
bcs.prop.2023
#pdf(paste0(path = "Lab_Data_TissueProcessing/output" ,"/BCS_proportion_2023.pdf"), height = 7, width = 13)
#print(bcs.prop.2023)
#dev.off()
#BCS_proportion_all+scale_fill_manual(values = c("#003C30","#01665E","#80CDC1","#C7EAE5","#F6E8C3","#DFC27D","#BF812D", "#8C510A","#543005" ))
#scale_fill_manual(values = c("#003C30","#01665E","#80CDC1","#C7EAE5","azure1","slategray1","slategray3", "slategray4", "gray25" ))
f <- function(pal) brewer.pal(brewer.pal.info[pal, "maxcolors"], pal)
(cols <- f("YlGnBu"))
#Coding graphs for LIS Conference - only 2023 data (continued)
df_3_prop<- select(data_2023, -c("site"))
#Count and proportion table for 2023 bcs scores by month
bcs_table <- table(data_2023$condition_score, data_2023$month, data_2023$site)
names(dimnames(bcs_table)) <- c("Body Condition Score", "Month", "Site")
bcs_table <- addmargins(bcs_table)
bcs_table
sessionInfo()
Add a new chunk by clicking the Insert Chunk button on the toolbar or by pressing Ctrl+Alt+I.
When you save the notebook, an HTML file containing the code and output will be saved alongside it (click the Preview button or press Ctrl+Shift+K to preview the HTML file).
The preview shows you a rendered HTML copy of the contents of the editor. Consequently, unlike Knit, Preview does not run any R code chunks. Instead, the output of the chunk when it was last run in the editor is displayed.